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ABSTRACT 



We present Gemini-S/GMOS-IFU optical spectroscopy of four regions near the centre of the nearby (3.8 Mpc) dwarf starburst galaxy 
NGC 5253. This galaxy is famous for hosting a radio supernebula containing two deeply embedded massive super star clusters, 
surrounded by a region of enhanced nitrogen abundance that has been linked to the presence of WR stars. We detected 1 1 distinct 
sources of red WR bump (Civ) emission over a 20'' (~350 pc) area, each consistent with the presence of ~1 WCE-type star. WC 
stars are not found coincident with the supernebula, although WN stars have previously been detected here. We performed a multi- 
component decomposition of the Ha line across all four fields and mapped the kinematics of the narrow and broad (FWHM=100- 
250 kms~0 components. These maps paint a picture of localised gas flows, as part of multiple overlapping bubbles and filaments 
driven by the star clusters throughout the starburst. We confirm the presence of a strong Ha velocity gradient over ~4V5 (~80 pc) 
coincident with the region of N/0 enhancement, and high gas density known from previous study, and interpret this as an accelerating 
ionized gas outflow from the supernebula clusters. We measure the ionized gas abundances in a number of regions in the outer IFU 
positions and combine these with measurements from the literature to assess the radial abundance distribution. We find that the 0/H 
and N/H profiles are consistent with being flat. Only the central 50 pc exhibits the well-known N/0 enhancement, and we propose 
that the unusually high densities/pressures in the supernebula region have acted to impede the escape of metal-enriched hot winds 
from the star clusters and allow them to mix with the cooler phases, thus allowing these freshly processed chemicals to be seen in the 
optical. 
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1. Introduction 

Local starbursts are the present-day manifestations of the intense 
star formation thought to have been instrumental in building 
galaxies in the early Universe. Starbursts occur in dense, gas- 
rich environments where the structure and kinematics of the ion- 
ized medium reflect the intense energy input from massive stars 
via stellar winds and supernova explosions. These winds can 
disrupt the surrounding gas, drive material out from the galaxy, 
and trigger (or quench) further star formation. 

NGC 5253 is an excellent nearby (distance of 3.8 Mpc; 
|Sakai et ar||2004| ) example of a starburst. This Blue Compact 
Dwarf (BCD) has ongoing active star formation at a rate o f 
-0.2 Mo yr"^ fealzetti et a l.|2007[|L6pez-Sanchez et a l.'2012'), 
and cont ains a large number of compact young (-1-12 Myr) star 
clusters ([Cresci et al.|[2QQ5t |Gorjian|p^96l [Harris et alT|[2QQ4| 
Vanzi et al.||2006 ) located within the central -250 pc starburst 



ters have been found in the outskirts of the galaxy underlining 
the fact that the current starburst is just the most recent star form- 
ing episode ( [Harbeck et al. 2012) . Lopez-Sanchez et al.|( [20T2 



concluded that the s tarburst was most likely triggered by the in- 
fall of an Hi cloud ( [Meier et al.|2002| ), although they could not 
rule out a contribution from ram pressure stripping as the galaxy 
moves th rough the dense intergalac tic medium of the Centaurus 
A group ( [Karachentsev et al.|2007 1. 

Ha images show a centrally concentrated region extending 
-30'' in the north-south direction, surrounded by a chaotic ion- 
ized structure with radial filaments, bubbles and loops ( |Calzetti| 
[et al. II 1999t [Marlowe et al. 1995; |Martin| | 1998 ). Dense clouds 
within these structures produce rapidly varying levels of extinc- 
tion. The morphology of the X-ray emission that traces the high- 
temperature gas follows the Ha distribution closely, and is con- 
sistent with a high energy-density outflow similar to that found 



zone. Optical studies have b een able to identif y only very few 
clusters older than 20 Myr ( [Harris e t^. 2004: Tre monti et aL] 
[2001 ) , supporting the extreme youth of this starburst. Recently 
three potentially massive (>10^ M©) and old (1-2 Gyr) star clus- 



* Based on observations with the Gemini-South telescope under pro- 
gramme GS-2008A-Q-25. 



in the dwarf starburst galaxy NGC 1569 ( [Strickland & Stevens 
[T9991 [Summers et al.|2004j . 

Within the region of brightest Ha emission are embedded 
two very massive (combined mass -1-2x10^ M©), young (age 
<4 Myr, but possibly <6 Myr; Martm-Her nandez et al. 2005]), 



compact super star clusters (SSCs) separated by -QV4 ( [Alonso- 



Herrero et al. ,2004^ (Vanzi & Sauvage,2004] ). These clusters ex- 
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ist behind Ay =7-1 7 magnitudes o f extinction ( [Alonso-Herrero 
|et al.|2QQ4{ |Vanzi & Sauvag e|2004| ), and are connected with (and 
possibly co-spatial with) a very strong mid-IR (Gorj ian et al. 
20011) and radio emission source ( [Turner &^B eck 2004 ; T urner 



et al.|200 ). Hereafter we will refer to this central giant Hn re- 
gion as the "supernebula" (following [Turner et al.|2000| . 

NGC 5253 is classified as a low-metallicity system (-0.2- 
0.3 Z©; Kobulnicky et al. 1999 ) with a fairly uniform oxygen 
distribution. Famously, however, a factor of 2-3 nitrogen en- 
hancement in the central regions has been found in a num- 



ber of targeted studies ( Kobulnic ky et al.|1997[ [Lopez-Sanchez 



[et al."2007i Wals h & Roy |, 1989 ). Evidence for the presence of 
nitrogen-rich WN- and carbon-rich WC-type Wolf-Rayet stars 



( Schaerer et al.|1997] ) in the same N-enriched regions led to the 
suggestion that these WR stars could be the cause of this en- 
hancement (|^bulmc^eni^ 1989). Indeed, 
in general, WR galaxies (those with large WR populations) as 
observed with SDSS samples exhibit elevated N/0 ratios rela- 
tive to non-WR galaxies ( [Brinchmann et al.|2008[ ). 

" Monreal-Ibero et aL] ( 2010 2012[ ) presented the results of 



the first spatially-resolved study of the ionized gas kinemat 
ics and abundances in the central regions of NGC 5253 with 
VLT/FLAMES integral field unit (IFU) observations. They 
found that, while the abundances of O, Ne and Ar remain con- 
stant within <0.1 dex over the mapped area, there are in fact 
two regions of elevated N/0: one strong enhancement associ- 
ated with the supernebula, and one weaker one associated with 
two older (~10 Myr) clusters towards the south-west. The fact 
that signatures of WR stars are not found at the location of this 
second enhancement led to the conclusion that there must be a 
more complex relationship between the presence of WR stars 
and N/0 excess than previously thought. Interestingly, [James 



et al. ( 2012b ) find high N/0 ratios in Haro 1 1 in a region with an 
integrated age just older than the WR phase (~7 Myr), but nor- 
mal N/0 ratios in two regions with large WR star populations, 
implying that the mixing time- scales are fairly long. 

In this paper we present optical IFU observations of four re- 
gions within the NGC 5253 starburst, including the supernebula 
region and two additional fields 150-250 pc away. We begin by 
presenting the distribution of carbon-rich Wolf-Rayet (WR) stars 
(WC-type stars; via the -5800 A red WR bump) in all four re- 
gions (Section [3]), since this has not been examined previously. 
We then present a detailed study of the kinematical properties of 
the ionized gas in all four regions (Section |4]) and the nitrogren 
and oxygen distribution in the two outer regions (Section|5]), and 
in comparison to measurements from the literature. A detailed 
study of the Hei abundance distribution within the central re- 
gions will be presented in a forthcoming paper (Monreal-Ibero 
et al. in prep.). Finally we discuss the state of the ISM across the 
whole starburst region, including the correspondence between 
the kinematics, chemical composition and presence of WR stars 
in relation to the enigmatic supernebula. 

2. Observations 

In 2008 queue-mode observations using the Gemini-South 
Multi-Object Spectrograph (CMOS) Integral Field Unit (IFU; 
Allington-Smith et"ar]|2002[ ) in one- slit mode were obtained of 



four regions near the centre of NGC 5253 (programme ID: GS- 
2008A-Q-25, PI: M. Westmoquette), in 0.5-0.7 arcsec seeing 
(-9-12 pc at the distance of NGC 5253) and photometric con- 
ditions. A nearby bright star was used to provide guiding and 



tip-tilt corrections using the GMOS on-instrument wave front 
sensor (OIWFS). 

The one- slit mode IFU gives a field-of-view (FoV) of 5 x 
3.5 arcsec (-90 x 60 pc) sampled by 500 contiguous hexagonal 
lenslets of 0V2 diameter. We took a set of dithered exposures 
at each position, off'setting by integer spaxel values (i.e. 0V2) 
and used the R831 grating to give a spectral coverage of 4750- 
6850 A at a dispersion of 0.34 Apix"^ An additional block of 
250 lenslets (covering 5 x 1.7 arcsec), off'set by 1 arcmin from 
the object field, provides a dedicated sky view. Positions 1, 2 
and 4 were tilted by 5° in order to avoid light from the galaxy 
from contaminating the sky fibres. Table [T] lists details of the 
observations for each position. 

The GMOS spectrograph is fed by optical fibres from the 
IFU which reformats the arrangement of the spectra for imaging 
by the detector. In order to remove the fibre-to-fibre through- 
put and detector pixel-to-pixel sensitivity diff'erences, and en- 
able the wavelength and flux calibration of the data, a number of 
bias frames, flat-fields, twilight flats, arc calibration frames, and 
observations of the photometric standard star Hiltner 600, were 
obtained as part of the standard baseline calibration plan. 

The IFU positions were chosen to cover selected areas of 
the disturbed ionized interstellar medium in the inner regions 
of NGC 5253, covering the nebulosity surrounding the central 
SSCs ( (Alonso-Herrero et al.[[2QQ4t [Turner et al.| |2000) and a 
further two regions to the south containing a number of ion- 
ized clumps. In Fig. [T] we show the position of the IFU fields 
on an HST/ACS-URC colour composite image (PID 10609, PL 
Vacca). 



2.1. Reduction 

The field-to-slit mapping for the GMOS IFU reformats the lay- 
out of the fibres of the sky to one long row containing blocks of 
100-150 object fibres interspersed by blocks of 50 sky fibres. In 
this way, all 750 spectra (500 object -\- 250 sky) can be recorded 
on the detector simultaneously. 

Basic data reduction was performed following the standard 
Gemini reduction pipeline (implemented in irai[j). Briefly, a 
trace of the position of each spectrum on the CCD was first pro- 
duced from the flat field. Then throughput correction functions 
and wavelength calibration solutions were created and applied to 
the science and standard star data, and the averaged sky spectrum 
(computed from the separate sky field) was subtracted, before the 
individual spectra were extracted using the flat-field trace. The 
spectra were then flux calibrated from observations of the stan- 
dard star. Although all our observations were taken under photo- 
metric conditions, that the standard star was not observed on the 
same night as our science observations means that the accuracy 
of the spectrophotometric flux calibration cannot be absolute. To 
test the accuracy of the relative flux calibration across the field, 
we re-reduced one science frame for IFU position 1 without sky 
subtraction (but with the standard fibre-to-fibre throughput cor- 
rection based on the twilight flat), and measured the flux of the 
prominent [Oi]/l5577 sky line across the whole field. Very little 
variation was found with no systematic patterns; the mean flux 
was measured to be 1.17 x 10"^^ erg s"^ cm"^ with a standard 
deviation of 8.80 x 10"^^ erg s"^ cm"^ (i.e. a -7.5% variation). 
This is at a level sufl&cient for our science goals. 



Position 1 was observed in exceptional seeing of 0V46. 



^ The Image Reduction and Analysis Facility (iraf) is distributed by 
the National Optical Astronomy Observatories which is operated by the 
Association of Universities for Research in Astronomy, Inc. under co- 
operative agreement with the National Science Foundation. 



Article number, page 2 of 



16 



M. S. Westmoquette et al.: Piecing together the puzzle of NGC 5253 




13^ 39^ 56.5^ 56.0^ 55.5^ 55.0^ 

Right Ascension 



Fig. 1. HST/ ACS-HRC colour composite (F330W, F435W, F550M, F814W, F658N) image of NGC 5253 s howing the location of the four 
GMOS-IFU positions. The dashed boxes indicate the location of the VLT/FL AMES -ARGUS IFU pointings from |Monreal-Ibero et aL| ( [20T0] ) and 
|Monreal-Ibero et ar]p012| ) for comparison. 

Table 1. Gemini/GMOS-IFU observation log 



Pos. 


Date 


Coordinates 


PA 


Exp. Time 


Seeing 


No. 




(J2000) 


o 


(s) 




1 


15/2/08 


13^39^56^0 -3r38'24V8 


355 


1500 X 4 


~0V5 


2 


1-2/3/08 


13^39^^55^.87 -3r38'24V7 


355 


1500 X 4 


~0V5 


3 


1/3/08, 9/7/08 


13^39^55^.72 -3r38'31V8 





1500 X 6 


~0V6 


4 


3/3/08, 9/7/08 


13^39^^55^.93 -3r38'36V9 


5 


1500 X 6 


0V5-0V7 



The data were converted into standard 'cube' format using 
GFCUBE with a spatial (over)sampling rate of OVl, and corrected 
for the effects of differential atmospheric correction (using an 
iRAF-based code written by J. R. Walsh based on an algorithm 
described in Walsh & Roy 1990 ), Finally, the individual dithered 
exposures for each position were combined using a custom pyraf 
script written by J. Turner (Gemini). 

In order to determine an accurate measurement of the in- 
strumental contribution to the line broadening, we fitted single 
Gaussians to isolated arc lines on a wavelength calibrated arc 
exposure, and found the resolution to be 1.2 A at Ha (55 km s"^) 
and 1.3 A at HjS (80 kms"^). The fit residuals shows that the 
instrumental line profile is extremely close to a Gaussian. 

Eig. |2] shows an example reduced and flux calibrated spec- 
trum extracted from a 13-spaxel region covering a bright ionized 



knot in lEU position 3 (which we define as region 1 in Section[5]). 
A plethora of emission lines are detectable from neutral to dou- 
bly ionized species, and in this spectrum the red Wolf-Rayet 
bump is also clearly s een. Line identificati ons were made using 
the lists compiled by [Tsamis et al.| ( |2003| ) and |L6pez-Sanchez| 
[etaLl ( |2QQ7| ). 



3. The carbon-rich Wolf-Rayet star population 

Because of the link posited in previous studies between WR 
stars and nitrogen enhancements, it is of particular interest to 
understand the population and distribution of WR stars within 
NGC 5253. Here we discuss our results on the WC star popula- 
tion, together with other WR star detections from the literature. 
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5000 5500 6000 6500 

Wavelength (A) 



Fig. 2. Example spectrum summed over 13 spaxels covering a bright ionized knot in IFU position 3 (defined as region 1; see Section |5]). All 
identified emission lines are labelled, including the broad C iv /15808 feature originating in the atmospheres of WC-type Wolf-Rayet stars. 



which together provide a view of the WC and WN star popula- 
tion within the central -350 pc of the galaxy. 

As exemplified in Fig. [2j we detect strong WR-star signa- 
tures in some spectra in the form of the red WR bump. This 
feature is a blend of the broad C iii /15698 and C iv /15808 emis- 
sion lines originating in the atmospheres of carbon-rich WC-type 
WR stars. At low metallicities such as in NGC 5253, early-type 
WC stars (WCE-types) are expected to be the dominant sub-type 
( |Crowther|20Q8| ). 

We mapped out the strength and location of the red WR 
bump in our four IFU fields. First we measured the contin- 
uum level in two windows either side of the WR feature and 
subtracted the average. We then collapsed the datacubes over 
the wavelength range 5765-5875 A to produce integrated maps 
of the red WR bump for all four positions. These are shown 
in Fig. [3] The size of each WR emission region is consis- 
tent with being a point source (i.e. the size of the seeing disk, 
-0.5-0.7 arcsecs). To estimate the number of WR stars in each 
source, we used the relation given by ( [Lopez-Sanchez & Este-| 
|ban|pnTO\ for the Civ ^5808 luminosity of a WCE star for a 
given metallicity (their equation 8). At the metallicity and dis- 
tance of NGC 5253, this relation predicts a WCE star luminosity 
of 2.3 X 10^^ erg s"^; given the uncertainties in both our mea- 
sured fluxes and the WCE star luminosity relation, each source 
is consistent with ~l-few WCE stars (see Table |4]k 

For a typical cluster mass of 5 x 10^ M© ( Har ris et al.|2004 ) 
and age of 4 Myr (Table [2]), the number of WR stars predicted by 
Starburst99 ( (Leitherer et al.|1999| ) is ~2. Firstly, this is in good 



agreement with our finding of ~1 WR star per source. Secondly, 
this very low number means that the presence of WR stars in any 
given cluster is dominated by stochasticity. Not every cluster 
of the right age and mass will contain a WR star (i.e. the IMF 
is not fully sampled in clusters of this mass). Indeed, simple 
Monte Carlo sampling of a Chabrier IMF over a total mass of 
5 X 10^ Mq, shows that only -1-10 stars >30 M© are expected. 

We do not detect a red WR bump at the location of the su- 
pernebula (the peak of emission in the optical; cluster H04-l|^, 
but find a number of WC stars within y including one coinci- 
dent with cluster H04-21, and a concentration at the UV emis- 
sion peak (cluster H04-4). We also detect four more WR emis- 
sion sources much further from the supernebula in GMOS posi- 
tions 3 and 4. 

[Schaerer et aL] ( [T997| detected W-R features (from both WN- 
and WC-type stars) at both near the peak of emission in the op- 
tical and the UV, equivalent to ~10 WC stars and a few 10s of 
WN stars in each of the two regions. Due to a pointing error 
during their observations they did not in fact cover the optical 
and UV flux peaks (clusters H04-1 and 4), but positione d their 
slit ~1V5 to the west (D. Schaerer, priv. comm. 2007). Sidoli 
( |2010| ), however, did obtain optical long- slit spectroscopy of the 
optical and UV peaks with VLT/UVES, and detected 5 WN and 



^ Hereafter we refer to the clusters identified in the study of Harris et al. 
(^004 ) with the notation H04-X, where X is their cluster ID. We note 
that their catalogue is by no means complete. Their study was based on 
WFPC2 data; more recent ACS-HRC imaging resolves very many more 
clusters (e.g. Fig.[3|. 
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1 WC star in a region covering the supernebula (cluster 

H04-1), and 10 WN and 2 W C stars at the UV peak (cluster 
H04-4). |L6pez-Sanchez et al.| 



(,2007 ) found ~1 WN star each 



in their regions A, B (a few arcsecs west and north-west of the 
supernebula), and ~10 WN stars in region C (equivalent to clus- 
ter H04-4). Accounting for the different distances and WR sub- 
type luminosities adopted by Lopez- Sanchez et al.| ( |2007] ), these 
values are in reasonable agreement with those found by ISidoli 
2010| ) in the corresponding areas. |Monreal-Ibero et a l. (2010) 



identified seven distinct regions containing blue-bump emission 
from WN-type stars in their IFU data covering the nuclear re- 
gions. The location of their detections are consistent with those 
found by the above authors. 

Fig. |4] shows a zoom-in of the central regions of NGC 5253 
summarising the IFU data on the location o f the WC star popula- 
tion from this study and the WN stars from| Monreal-Ibero et al.| 
( [2010| ). In only two locations are there coincident detections of 



WC and WN stars (regions WR3 and H ii 2 of Monreal-Ibero 
et al.pOlO] ), although the spatial overlap of our datasets is far 



from complete. Of the four regions in which nebular He ii (i.e. 
narrow emission, as opposed to the very broad emis sion orig- 
inating in the WR stellar atmosphere) is detected by [Monreal-] 



Ibero et al. (2010), two are covered by our GMOS data, and both 



contain a WC star (regions He ii 1 and He ii 4 of |Monreal-Ibero 



|et al.pOl O^; see Fig. |3]). Only certain classes of WR stars, and 
possibly certain O stars, are predicted to be hot enough to ion- 
ize nebular He ii ( [Crowther||2007 [ [Kudrit zki 2002 ), whose ion- 
ization potential is hv = 54.4 eV, and so far only one WC star 
has been associated with a He ii nebula (WC star MC45 in neb- 
ula BCLMP38b in M33; |Kehrig et al.|2011| ). Our spectroscopy 
shows that both our He n-WC regions have higher than average 
[O in]/H/5 line ratios, consistent with the presence of a hard ra- 
diation field. Neither, however, show convincing evidence for 
the presence of shocks in the traditional EFT diagrams (Bald win] 
let al.|1981| not shown). 

Almost all previous studies of the WR population in 
NGC 5 253 have focussed on the central regions. However, |Sidoli| 
([2010') examined a region ~5'' to the nor th of the supernebula, 
and found evidence for 1 WN star, and [Lopez-Sanchez et al.| 
( [2007| found evidence for a further WN star in their region D 
(coincident with cluster H04-2, -10'' south of the supernebula 
and ~5'' to the east of GMOS IFU position 3). We detect red- 
bump WR emission consistent with the presence of single WC 
stars in 4 locations in GMOS positions 3 and 4. In position 3, 
one of these WC stars is coincident with H04 -17, a known clus - 
ter of mass 2-3x10^ M© and age <4 Myr (Harris et al.||2004t . 
The remaining 3 sources are not coincident with any previously 
identified clusters. Table [2] summarises the detections of WN 
and WC stars in the regions defined in both this study and that 
( [2QTq1 ) 



of Monreal-Ibero et al, 



where we have ordered them by 
decreasing EW(HyS) (or increasing age using the predictions of 
Starburst99). 

For all stars that pass through a WR phase, it is thought that 
the WN stage precedes the WC stage ( Crowther^2007j . Thus, 
one would expect WN stars to be present predominantly in the 
younger regions, and WC stars in the older regions. Disregard- 
ing regions 2 and 6 (where we do not detect a red WR bump, and 
that were not covered by the observations of [Monreal Ibero et^L] 
[2010 ), there does appear to be a transition at ~ 3.4-4.4 Myr, be- 
low which there is a higher probability of finding WN stars, and 
above which there is a higher probability of finding WC star. 
Clear exceptions to this are regions H ii 2 (very young but ex- 
hibits evidence for both WR types) and WR5 (which, according 
to its UJ3 flux, is older but has WN star signatures). In neither of 




Fig. 4. Zoom in to the central region of NGC 5253 showing all the 
available IFU data on the location of the WR star population. Within 
GMOS IFU positions 1, 2 and 3, the red WR bump flux is represented 
with white contours. The size of the seeing disk for positions 1 and 
2 (0^5) is shown in the top left. The large whit e dashed rectangle 
shows the VLT/FLAMES -ARGUS IFU position of Monreal-Ib ero et al.| 
( |2010| ), and the labelled cyan and magenta rectangles show the regions 
in which evidence for a blue WR bump and He ii emission was found, 
respectively. 



these locations did [Harris et aL] p004[ ) identify a cluster, but as 
we have previously mentioned, their catalogue is far from com- 
plete. Furthermore, Mon real-Ibero et al.[ ([2010) were only able 
to determine the location of any WR emission to a precision of 
~1V5, which is not sufficient to be sure which source (if identi- 
fied) is responsible. Projection eff'ects also add to the difficulties 
of disentangling which source is associated with which. 

In summary, we detect 11 distinct sources of CiV/15808 
emission throughout the NGC 5253 starburst, each consistent 
with the presence of ~1 WCE type star. Some are associated 
with known young clusters catalogued by H arris et al.[ ( [2004[ al- 
though this study is by no means complete), and two are coinci- 
dent with the location of WN stars detected by [Monreal-Ibero 



et al. 



( [2010[ ). The finding of WR stars spread out over 20 
(-350 pc) attests to the large area over which the starburst has 
occurred. 



4. Emission line kinematics 

Previous long-slit observations have shown that the ionized gas 
kinematics in the central regions of NGC 5253 are complex, 
with line profil es revealing multiple components and asynmiet- 
ric wings (e.g. [Lopez-Sanchez et al.[[2007[ [Martin & Kenni- 



cutt|1995t|Sidofi||2010||Vanzi et al.||2006| [Walsh & Roy||1989| ). 



Monreal-Ibero et al.[ ( [2010[ ) present a detailed analysis of the 
emission line kinematics in the central regions as covered by 
their VLT/FLAMES -ARGUS observations. They found two- 
three components were needed to fit the observed Ha line pro- 
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13^ 39^ 56.5^ 56.0^ 55.5^ 13^ 39^ 56.5^ 56.0^ 55.5^ 

Right Ascension Right Ascension 



Fig. 3. Left: Integrated maps of the C iv /15808 and C iii /15698 WC-star emission feature (the red WR bump) for each position on the same colour 
scale (as shown), overlaid on an //5r/ACS-HRC F550M image (ID: 10609, PI: Vacca). The emission flux in each knot (consistent with a point 
source the size of the seeing disk) is consistent with ~1 WCE-type star. The HRC F550M image is reproduced on the right with a number of the 
clusters identified by Harris et al. ( 2004 ), and discussed in the text, labelled. 

Table 2. Yip equivalent wi dths (EW), ages and presenc e of WR stars in regions from both this study (bold numbers; see Table [J] for details of the 
measurements) and that of |Monreal-Ibero et al.| ( |2010) . 



Region 


SN^ 


Hii2 


WR2 


WRl 


WR3 


WR4 




5 


6 


1, He II 4 


WR5 


2 


4 


3 


Hen 1 


Other literature 


H04-1, 








H04-4, 




H04-21 






H04-17 












names^ 


LS-B 








LS-C 






















EW(Hy6) (A) 


250-350 


320 


180 


130 


130 


130 


110 


105 


100 


95 


85 


60 


34 


14 


9 


Age (Myr)'^ 


2.3-2.7 


2.4 


3.2 


3.4 


3.4 


3.4 


4.4 


4.5 


4.6 


4.6 


4.7 


4.8 


5.4 


6.3 


7.5 


WN'^ 


/ 






y 


y 


y 


X 


? 


? 


X 


y 


7 


7 


7 


X 




X 


y 


? 


? 


y 


? 


y 


y 


X 


y 


7 


X 


y 


y 


y 
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file, and attributed them (at least near the central SSC) to cluster 
winds and multiple expanding shells. 

As mentioned above, |Sidoli| ( [2Q1Q| ) presented VLT/UVES 
long- slit echelle spectroscopy of the optical and UV peaks (knots 
1 and 2 from the study of |Monreal-Ibero et al.|2010| . They also 
found three Gaussian components were required to satisfactorily 
fit the emission line profiles: a narrow (FWHM=30-40 kms"^) 
component, a broad (FWHM=100-150 kms"\ blueshifted by 
10-30 km s"^ with respect to CI) and a third, fainter component 
intermediate in width and radial velocity. With a spectral res- 
olution of R~28 000 (~1 1 km s"^ at Ha), the fact that they still 
needed a broad component suggests that it is not simply a convo- 
lution of many narrower, unresolved components (c.f. |Melnick| 
fet al.,1999, for 30Dor). 

4.1. Decomposing the line profiles 

We fitted multiple Gaussian profile models to the Ha emission 
line using an idl;^^ fitting package called pan (Peak ANalysis; 
Ipimeo 2005 ), to recover kinematic information about the ion- 
ized gas in the four IFU positions. The high spectral resolu- 
tion and signal-to-noise (S/N) of our data, together with the fact 



that the instrumental line profile is so close to a Gaussian (Sec- 
tion|2.1[), have allowed us to quantify the Ha line profile shape to 



a high degree of accuracy. We can therefore be confident of the 
identification of faint components and/or those with very small 
(<10 kms"^) velocity separations. 

To account for the variable number of Gaussian components 
in each line profile, we fitted a single, double, and triple Gaussian 
component initial guess to each. Line fluxes were constrained to 
be positive and widths to be greater than the instrumental con- 
tribution to guard against spurious results. Multi-component fits 
were run several times with diflferent initial guess configurations 
(widths and wavelengths) and the fit with the lowest value 
was kept. However, we note that the minimisation routine 
employed by pan is very robust with respect to the initial guess 
parameters. 

To determine how many Gaussian components best fit an ob- 
served profile (one, two or three), we used a likelihood ratio test 
to determine whether a fit of ^ -I- 1 components was more appro- 
priate to an ^-component fit. This test says that if the ratio of 
the x^ statistics for the two fits falls above a certain threshold, 
then the fits are considered statistically distinguishable, and the 
one with the lower x'^ can be selected. Here we determine the 
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threshold ratio by visual inspection of a range of spectra and fits. 
This method is a more generalised form of the formal F-test that 
we have used in previous work, e.g. Westmoquet te et al.| pQQ9, 
2Q11| ). We are aware of work that cautions against the use of the 
F-test to test for the presence of a line (or additional line com- 
ponent; e.g. |Protassov et al.|[2002] ), but given the absence of a 



statistically correct alternative that could be applied sensibly to 
the volume of data presented here, we have chosen to opt for this 
generalised;^^ ratio test approach. 

This test, however, only tells us which of the fits (single, 
double or triple component) is most appropriate for the corre- 
sponding line profile. Experience has taught us that we need to 
apply a number of additional, physically motivated tests to filter 
out well-fit but physically improbable results. We specified that 
the measured FWHM had to be greater than the associated er- 
ror on the FWHM result (a common symptom of a bad fit), the 
fluxes of all components should be >0, and we rejected any fits 
where = (another symptom of a bad fit). 

The final step of the fitting procedure was to assign each fit 
component to a particular map in such a way as to limit the con- 
fusion that might arise during analysis of the results, such as 
discontinuous spatial regions arising from incorrect component 
assignments. In positions 1, 2 and 3, we found that the maps 
were most coherent when we assigned the broadest component 
to C2, and after that, the brightest as CI and the faintest to C3. 
In position 4, we simply assigned the brightest component to CI 
(regardless of width), then the faintest component to C3. This 
approach helps to limit the confusion that might otherwise arise 
during analysis of the results where discontinuous spatial regions 
might arise from incorrect component assignments. 

Experience has also shown us that the errors that pan quotes 
on its fit results are an underestimate of the true uncertainties, 
particularly once uncertainties on our post-fitting tests and filters 
are taken into account. We can therefore estimate more real- 
istic errors through a visual re-inspection of the profile-hfit af- 
ter knowing which one was selected by our tests, and by taking 
into account the S/N of the spectrum. We estimate that errors in 
FWHM range between 0.25-2 kms"^ for CI, 4-20 kms"^ for 
C2 and 15-25 kms"^ for C3, whereas those for radial velocities 
range between 0.1-3 kms"^ for CI, 2-5 kms"^ for C2 and 10- 
30 kms"^ for C3. A 5-10 percent increase in the CI and C2 
uncertainties may occur with the addition of a third component. 

Figs. |5j [6] and [7] show the resulting FWHM and radial ve- 
locity maps for positions 1, 2, 3 and 4, respectively. In these 
figures we also include a number of example Ha line profiles 
and best-fitting Gaussian models from certain spaxels. These ex- 
amples were selected to illustrate the discussion below, but also 
demonstrate the high quality of the spectra and the accuracy of 
the line-fitting. 

4.1 .1 . Positions 1 and 2 



As mentioned above, | Monreal-Ibero et aL] ( [2010| ) present a de- 
tailed analysis of the emission line kinematics in the central re- 
gions. However, the higher S/N and spatial resolution of our ob- 
servations permits us to examine the emission line profiles and 
kinematical morphologies in the central regions in more detail 
(positions 1 and 2 were taken in seeing of -0.5'')- We present 
our results in the following. 

IMonreal-Ibero et al.| ( |20T0l ) and [Sidolil (20W) find broad 
line components in the supernebula region with widths of 
FWHM~80-150 kms"^. In our analysis, we have been able to 
identify components with FWHM>250 kms"^ in a handful of 



spaxels. These very broad components are faint (~5 percent of 
the total flux); an example is shown in Fig.jSJ). In these spax- 
els, the need for this additional component is very clear from 
the residuals and of the fit. In most of the field, however, the 
broad component line widths are consistent with those found in 
the aforementioned studies (FWHM~ 100-1 50 kms"\ compris- 
ing 10-50 percent of the total line flux). 

|Monreal-Ibero et aL] ( |2010| ) identified a velocity gradient in 
the broad component extending over ~4V7 (-85 pc) centred on 
the supernebula and oriented NW-SE. This gradient is clearly 
seen in our C2 velocity map (Fig. |5] bottom centre panel), and 
there are clear hints of this gradient also in the secondary nar- 
row component (C3). Our IFU map shows that this gradient is 
present across the entire C2 map, with an appearance similar to 
a rotation field. However, there are distinct peaks on the red and 
blue side, hinting that there may be other explanations. 

To investigate this further, we defined a 0V7-wide pseudoslit 
crossing the supernebula, intersecting the two peaks, and ori- 
ented perpendicular to the line of nodes tracing zero velocity 
(i.e. along the maximum velocity gradient, as shown in Fig. [5]). 
From this we extracted the radial velocities of the three line com- 
ponents to plot in Fig. [8] From this figure, it is clear that whereas 
CI, the dominant narrow component, remains fairly static across 
the whole region with velocities varying by only ~30 kms"\ 
C2 and C3 trace a linear gradient in velocity from — 50 to 
~-h50 kms"^ within the central ~4'' (70 pc). Beyond this, C2 
extends to even higher velocities, but is not followed by C3. All 
three components cross at zero velocity within ~0V5 of the loca- 
tion of the supernebula. Rodr iguez-Rico et al.| (|2007) presented 
observations of the radio recombination line H53Qr over a similar 
region, and found a -50 kms"^ gradient over ~3V5 in the same 
NW-SE direction. That this is in very good agreement with our 
Ha measurements indicates that the ionized material throughout 
the gas column in this region all behaves similarly. 

Monreal-Ibero et alT] ( |2010| ) interpreted this velocity gradient 



(seen by them only in the broad component) as tracing an out- 
flow from the SSCs in the supernebula, and our data here also 
support this conclusion. The fact that the velocity gradient exists 
in the broad line component is difficult to reconcile with rotation 
since such a turbulent rotating structure could not maintain co- 
herence for long. We therefore conclude that C2 and C3 trace 
an accelerating ionized gas outflow (although we cannot rule out 
the possibility that the outflow is superimposed on a more gen- 
eral rotation-like bulk motion existing separately to the majority 
of the ionized gas traced by component CI). The fact that the su- 
pernebula is located very close to the zero-point of the velocity 
gradient implies that the clusters embedded there are the driving 
sources of the outflow. 

Proceeding with the outflow interpretation, its inclination, /, 
is still unknown; assuming /=45° the deprojected maximum out- 
flow speed is ~70 kms"^ This velocity is roughly equivalent 
to the line width of C3 (FWHM~30-100 km s"^), but is signifi- 
cantly less than the width of C2 (FWHM~100-250 km s"^). This 
suggests that the turbulent energy in the outflow exceeds its bulk 
flow energy. 

4.1.2. Position 3 

IFU position 3 covers a region -6-10'' to the south-west of the 
supernebula containing two Hii knots, the eastern of which is 
associated with a known cluster of mass 2-3x10^ M© and age 
<4 Myr (H04-17). A number of other young star clusters not 
catalogued by |Harris et al.| ( [20Q4| ) are also located in the field. 
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Fig. 5. Pos 1 and 2 Ha flux (top panels), FWHM (middle panels) maps in units of km s~^ (corrected for instrumental broadening), and radial 
velocity maps (bottom panels) in units of km s~\ relative to Usys = -390 km s~^ ( Schwartz & Martin 2004') for the three identified line components. 
Contours show the summed Ha intensity (the peak at coordinates (3,3) is the supemebula). Example Ha line profiles and their best multi- 
component fits are shown for select spaxels as discussed in the text. The parallel dashed lines show the position and width of the pseudoslit we 
defined to extract the velocity information for plotting in Fig. [8] The dot-dashed line in the C2 radial velocity map (bottom-centre) shows the zero 
velocity (i.e. Usys) contour. 
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Fig. 8. Position- velocity diagram along the pseudoslit crossing the su- 
pernebula, and oriented along the maximum velocity gradient, as shown 
in Fig.[6| Velocities are shown relative to Usys = -390 km s~^ ISchwartz' 
[& Martin 2004 ). The vertical dashed lines indicate the region within 
which both C2 and C3 follow the same steep gradient. We interpret this 
as evidence of an outflow from the supernebula clusters. 



C2 (the broad component) is identified across the entire field, 
and exhibits line widths extending to FWHM>200 kms"^ The 
broadest line components are generally also blueshifted (exam- 
ple line profile (b) in Fig. [6]). Around the western Ha knot, the 
C2 component is fainter, narrower (-100 km s"^), and redshifted 
compared to the surrounding material (-30 km s"^); see example 
line profile (a). The corresponding CI component is blueshifted, 
suggesting the gas is expanding or outflowing from this Ha knot. 
To the south of the IFU field, we see regions of clear line splitting 
where the two components are more equally matched in terms of 
flux and width (example c). 

This position partially overlaps with the IFU field-of-view of 
|Monreal-Ibero et aL] ( |201Q| ), and here our CI radial velocities are 
in good agreement with their measurements. They find similar 
regions around their knot 3 in the south-west of their field where 
two distinct narrow components are seen. They conclude that the 
environment here (clusters H04-3 and H04-5; see Fig. [3]) must be 
more evolved than around the supernebula, and that these older 
clusters have cleared out their gas, leaving broken shells ionized 
by the remaining hot stars. 

4.1.3. Position 4 

Position 4 covers an interesting arc-shaped feature seen in Ha 
~12V5 (220 pc) to the south of the supernebula. There are a 
handful of young clusters in this area that appear to be located 
mostly along the rim of the Ha arc, including one identified by 
[Harris et al.| ( |2QQ4l H04-27). 
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Fig. 6. Pos 3 Ha flux (top panels), FWHM (middle panels) maps in units of km s ^ (corrected for instrumental broadening), and radial velocity 



maps (bottom panels) in units of km s ^ , relative to Usys = -390 km s ^ ( Schwartz & Martin 2004|) for the three identified line components. Contours 



show the summed Ha intensity. Example Ha line profiles and their best multi-component tits are shown for select spaxels as discussed in the text. 



The very broad Ha component seen in much of position 3 
is not seen here. Line widths mostly remain below 100 kms"^ 
across the field, with only a few regions where line widths extend 
up to -150 kms"^ The Ha line shape at the peak of emission 
is consistent with a single Gaussian. Across the rest of the arc, 
the line shape is characterised by a double-Gaussian (examples 
(b) and (c)) with velocity diflTerences between the components 
of 20-30 kms"^ In the north of the field, C2 is fainter and is 
blueshifted with respect to CI by 20-30 km s"^ These split line 
components indicate low- velocity flows of gas, likely to be shells 
or bubbles driven by the nearby star clusters. The bright Ha knot 
located at the northern edge of the field does not appear distinct 
in the kinematics maps. 

4.1 .4. Summary of kinematics 



Together with the results of |Monreal-Ibero et al.| p010), our re- 
sults paint a picture of the overall gas dynamics in NGC 5253 



consisting of localised gas flows, as part of multiple overlap- 
ping bubbles and filaments driven by the star clusters through- 
out the starburst. Close to the supernebula, we have been able 
to identify faint (~5 percent of the total flux) components with 
FWHM>250 kms"\ although over most of positions 1 and 2, 
the broad component line widths are consistent with those found 
in the aforementioned studies (FWHM- 100-150 kms"\ com- 



prising 10-50 percent of the total line flux |Monreal-Ibero et al 
|2QTQl|Sidoli|2QTQl ). 

We confirm the presence of a strong velocity gradient in the 
central region seen by |Monreal-Ibero et aL] ( |2010| ), extending 
over ~4V5 (-80 pc) oriented NW-SE direction. The p-v diagram 
extracted along the maximum gradient shows that both the broad 
and secondary narrow Ha components (C2 and C3) both trace 
this steep gradient, which crosses zero vel ocity at the loc ation of 
the supernebula. In line with Monreal- Ibero et al.| ( |2010| ), we in- 
terpret this as an ionized gas outflow from the supernebula clus- 
ters. Although the clusters are known to be deeply embedded 
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Fig. 7. Pos 4 Ha flux (top panels), FWHM (middle panels) maps in units of km s~^ (corrected for instrumental broadening), and radial velocity 
maps (bottom panels) in units of km s~\ relative to fsys for the three identified line components. Contours show the summed YLa intensity. Example 
Ha line profiles and their best multi-component fits are shown for select spaxels as discussed in the text. 



(Ay=7-17 mag; " Alonso-Herrero et al."2004i "Vanzi & Sauvage" 
|2QQ4 ) , this implies that the outflow has already carved channels 
in the dense gas distribution through which the ionized mate- 
rial is escaping. Furthermore, the fact that we see a continuous 
velocity gradient in both the blueshifted and the redshifted side 
implies that the extinction in the supernebula must fall off' rapidly 
such that both the front- and back- sides of the outffow are visible 
in the optical. 

In our IFU position 3, 6-1 C to the south, the broad compo- 
nent contributes <15 percent, and by position 4 is not detected 
at all. This suggests that its presence is connected to the con- 
centration of young clusters in the supernebula region, and most 
likely to the aforementioned outffow. However, since we de- 
tect broad emission >150 pc from this cluster it is unlikely to 
be solely due to its winds. It is possible that other young clus- 
ters are blowing winds that produce this broad line. Since the 
broad lines in the central regions follow the bulk motions of the 
more quiescent gas, whatever causes this turbulence must be a 
local phenomenon. [Westmoquette et aL] ( |20Q7| ) detected a very 
similar broad line component to the ionized gas emission lines 



in NGC 1569. Here we found the broadest emission correlated 
with the location of gas clumps in the ISM, and hence concluded 
that it traces turbulent mixing layers (TMLs; [Binette et al.|2QQ9t 
|Slavin et"aL]|1993| ) on the clump surfaces set up by t he impact 
of stellar/cluster winds. In [Westmoquette et al.| (|2Q10|) we tested 
this hypothesis by observing a gas clump near to the Galactic 
OB association Pismis24 (NGC 6357) at high spatial and spec- 
tral resolution. That we found broad emission originating from 
the clump surface in projection conffrmed the idea that TMLs 
can explain this kind of broad component emission. 

5. Chemical abundances 

Gas phase abundances in NGC 5253 have been of partic ular in- 
terest since the connection between N- enriched regions (|Camp-| 



bell et al.| 19861 [Kobulnicky et al.| 1 997 [ [Walsh & Roy|1989| ) and 

the presence of WR stars (WN-types in particular) was made 
d Schaerer et al.^1997) . However, it is only recently that an unbi- 
assed mapping of the abundances could be made using integral 
ffeld spectroscopy ( [Monreal-Ibero et al.|201Q|[2Q12| ). Here, they 
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represents a significant advantage over long-slit/fixed aperture 
studies, since positioning uncertainties and slit/aperture losses 
are avoided, and the measurements can be viewed in context 
since the surrounding environment is known. The location of 
our eight regions is shown in Fig.|9] 

Region 1 (~7V4 or 135 pc from the supemebula) corresponds 
to cluster H04-17 and to region He ii 4 from ( Monre al-Ibero et aL] 
2010). The full spectrum from this region is shown in Fig. [2j 



Region 2 corresponds to an Ha knot 1V8 to the west; region 3 
contains WR star signatures and likely corresponds to a cluster 



visible on the HST image that was not catalogued by Harri s et al. 
2004[ see Figs. |3] and |4]); region 4 exhibits WR signatures, but 



no strong Ha and there is no obvious cluster at this location in 
the HST imaging; region 5 correspon ds to an Ha knot and most 
likely to a cluster not catalogued by [Harris et"aL] p004| ); and 
region 6 corresponds to the Ha knot associated with an arc or 
bubble. We also defined two larger boxes covering regions of 
fainter, diff'use gas, d_3 and d_4, in order to sample the abun- 
dances away from the bright H ii knots. 




arcsecs 



0.2 0.4 0.6 0.8 1 

log(Flux) 

Fig. 9. IFU positions 3 and 4. Shown in grey scale is the integrated Ha 
flux. Red contours show the integrated red WR bump flux (dominated 
by C iV/15808 emission) originating from WC-type Wolf-Rayet stars. 
The black outlines show the 6 regions from which we extracted summed 
spectra for our abundance analysis. 



studied a -9x5 arcsec area encompassing the central regions. By 
mapping out the gas phase abundances at a resolution of ~0V8, 
they found two locations with enhanced N/0: one significant en- 
hancement associated with the supernebula region and another 
much smaller enhancement associated with two clusters ~7'' to 
the south-west (clusters H04-3 and H04-5). 



5.1. Definition of regions 

Our integral field observations provide the opportunity to in- 
vestigate the gas phase abundances much further from the su- 
pernebula. However, the ionized gas surface brightness falls 
off' rapidly with distance, and our observations are not of high 
enough S/N to derive spaxel-to-spaxel abundance maps as was 
done in Monreal-Ibero et al.|([2012 ). Thus we have used the Ha 
and WR ffux maps to define precise apertures from which we 
extract integrated spectra and perform our analysis. This already 



5.2. Calculation of abundances 

Abundances were derived using neat (Nebular Empirical Anal- 
ysis Tool; Wesson et al. 2012|j a code for calculating chemical 
abundances in photoionized nebulae. The code carries out a stan- 
dard analysis of lists of emission lines using long-established 
techniques to estimate the level of interstellar extinction, cal- 
culate representative temperatures and densities, compute ionic 
abundances from both collisionally excited lines and recombina- 
tion lines, and finally estimate total elemental abundances using 
an ionization correction scheme ( [Kingsburgh & Barlow||1994| . 
Additionally, neat uses a Monte Carlo technique to robustly 
propagate uncertainties from line ffux measurements and ioniza- 
tion correction factors (ICFs) through to the derived abundances. 
We note here that this more accurate treatment of the uncertain- 
ties produces reliable errors on the calculated abundances, and 
as such they may be larger than what can be found for equivalent 
measurements in the literature. 

We measured the integrated line ffuxes of all the major emis- 
sion lines in each of the extracted integrated spectra and list 
them in Table [3] We made no attempt to decompose each emis- 
sion line into individual components for this part of the analy- 
sis. Using these values as inputs into neat we derived the elec- 
tron temperature, Tq, the electron density, A^e, and ionic and 
elemental abundances for all six regions. The code was run 
with -10000 iterations, a number sufficient to establish a sta- 
ble uncertainty distribution (We sson et al.^2012 ). Tq's were de- 
rived from the temperature- sensitive [Nii]/15755//16584 line ra- 
tio, whilst A^e's were derived from the [S ii]/l6716//16731 line 
ratio. For Region 4, where [Nii]/15755 was undetected, we 
adopted the average of Regions 5 and 6, 9700 K. Since the 
[Oii]/13727 line was not within our wavelength coverage we 
adopt ed measuremen ts from the nearest corresponding regions 
from (Walsh & Roy| ([W89 ). For our regions 1-3 we used the 
va lue I([0 ii])/I(Hyg)=3 .083+0.034 (corresponding to 'region 2' 
of Walsh & Roy 1989), whilst for our regions 4-6 we adopted 
I([Oii])/I(HyS)=2.929±0.031 ('region 7' of Walsh & Roy 1989] ). 
As mentioned above, neat adopts the ionisation correction fac- 
tors (ICFs) listed by ^Kingsburgh & Barlow] ( |1994| ) to convert 
from ionic to elemental abundances, and propagates the statisti- 
cal uncertainties associated with the ICFs accordingly. Derived 



Full documentation, including source code and reference files, is 
available at http : //www . sc . eso . org/~rwesson/codes/neat/ 
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O, N, S, and He nebular abundances and their uncertainties are 
given in Table |4j 

Our derived O abundances are, within the uncertainties, con- 
sistent with a metallicity of -0.2 Z©, as has been found in other 
regions of the galaxy (for a discussion see below). One point lies 
significantly above the others (for Region 5), but has large un- 
certainties that derive from the low S/N of the [N ii]/l5755 line in 
this region. The N abundances suggest a slightly lower metallic- 
ity. The N abundance in these outer regions is significantly lower 
than that found in the super nebula region, which is enhanced to 
-0.4 Z© ( [Monreal-Ibero et al. 2012). O ur He/H values are con- 
sistent with what was found by Lop ez-Sanchez et aL] ( |2007 ). 

It is interesting to note how similar the spectra of the diff'use 
regions are to the emission knots/Hii regions. This implies that 
a significant fraction of the ionizing photons are leaking out of 
the gas clouds around the clusters to ionize the intra- [star]cluster 
medium. 



5.3. Radial abundance patterns 

We have compiled all the 0/H, N/H and N/0 abundance mea- 
surements for diff'erent regions found in the literature, and, to- 



gether with our new measurements, listed them in Table 5.4 



Fig. [T0| shows these values plotted versus distance from the su- 
pernebula (shown in Fig. [T] with a white cross). Note that one 
point at a particular radius is not necessarily associated with an- 
other point at the same/similar radius. Also included are the 



individual spaxel values from Monreal- Ibero et aL] ( |2012[ their 
fig. 8). As mentioned above, (within the uncertainties) the 0/H 
distribution is consistent with a metallicity of -0.2 Z©. Further- 
more, we can see that it does not show any systematic trend (in- 
creasing or decreasing) within the 250 pc area sampled. The 
N/H profile, however, shows the well-known enhancement of 
-0.5 dex in log(N/H) within the central 50 pc. The large scat- 
ter of points within the rise result from the fact that the N- 
enhanced region is not circular in shape ( [Monreal-Ibero et aL] 
2012[). At la r ger ra dii, both our measurements and those of 
Walsh & Roy (1989]^ are systematically lower than those mea- 
sured by [Monreal-Ibero et aL]p012| ), possibly due to the system- 
atic uncertainties introduced by the diff'erent temperature- and 
density- sensitive lines used. However, given these uncertain- 
ties, the plot is consistent with a flat N abundance beyond the 
supernebular region (i.e. R>50 pc). The N/0 distribution also 
reflects the N/H enhancement in the supernebula region, and is 
also consistent with being -flat beyond 50 pc. 

Over the years a number of low-mass dwarf galaxies and gi- 
ant Hii regions have been the subject of spectroscopic chemical 
abundance mapping (e.g.[Devost et a l. 1997; Gonzalez-Delgado[ 
et al.ll994|potov & Th uan'1999[ [Kobulnicky & Skillman|1996[ 



1997[ Russell & Dopita 1990), and more recently with IFUs 



(e.g. [James et al. 20101 |2012a[ ^ [2QQ9l [Lagos et al.[[2QT2l [Pgrez^ 



Montero et al.,2011. ). With very few exceptions (see below). 



these systems do not show significant internal chemical fluctua- 
tions or gradients, even though the expected yield of heavy ele- 
ments produced by massive stars in young starbursts is substan- 
tial. [Kobulnicky & Skillmanl ([1997 ) put forward two explana- 
tions: either ejecta from stellar winds and supernovae are trans- 
ported throughout the galaxy and mixed instantaneously and uni- 
formly, or freshly synthesised elements remain unmixed with the 



The uncertainties on the 



Walsh & Ro^ \m9\ measurements are 



small because they summed over much larger regions to get a high 
S/N and to explore possible similarities with more classical giant H ii 
regions. 



surrounding interstellar medium and reside in a hard-to-observe 
hot 10^ K phase (or a cold molecular phase). In support of the 
latter scenario, Martin et al.|(^2002 ) find evidence for a metal en- 
riched diff'use X-ray phase in the halo of the starbursting dwarf 
galaxy NGC 1569. Measurements like this, however, are ex- 
tremely hard to make with existing X-ray technology, and there 
are no equivalent X-ray measurements available for NGC 5253. 

NGC 5253 is one of the most well-known (and certainly 
nearest) examples amongst the rare cases of dwarf galaxies with 
significant chemical variations. These include Mrk996 (Tames 
et_al. 2009), UM448 (James et al. 2012a) and Haroll (Jan^ 
et al.|[2012b|), possibly ICI O ([Lopez-Sanchez et al. 2011) and 
IIZw70 (iKetiri g et al.|2008[), and those with radial oxygen abun- 
dance gradients ^Lee et al. [20061 [Werk al.[2010| ). Of these, it is 
only NGC 5253 and Mrk996 that show a genuine N/H and N/0 
enhancement (a high N/0 ratio by its self can result from the dilu- 
tion of the O abundance by an inffow of metal-poor gas). How- 
ever, although they both show these enrichments, the elevated 
N/H ratios in Mrk996 are only found in the broad emission-line 
component (FWHM>300 kms"^), and the associated gas densi- 
ties are even higher (10^-10^ cm"^). So what is it in NGC 5253 
that puts it in this extremely rare class of galaxies? 



5.4. What is the origin of the N/0 enhancement in the 
supernebula region? 

Putting together the available facts we know that: 



- The radio supernebula ( [Turner & Beck[[2004| ) contains two 
very massive (combined mass -1-2x10^ M©), young (ages 
<4 Myr, but see below), embedded super star clusters (SSCs) 
([Alonso-He rrero et al . 2004; Mart m-Hernandez et al.|2005[ 
"t^nzi & Sa uvage|20Q4| ). 

- N/0 is strongly enhanced in a region measuring -2''x4'' ex- 
tending in the NW-SE direction, and is not centred exactly 
on the supernebula cluster, but -0.5-1'' to the NW. 

- There does not appear to be any obvious relation between 
the N-enhancements and the ages of the clusters and the WR 
populations. WR stars of both varieties are found up to 10'' 
from the SN, but there is no corresponding N enhancement at 
these locations. Apart from the supernebula region, only one 
other, much weaker, N/0 enhancement was found around 
clusters H04-3 and 5 ( [Monreal-Ibero et al.[[2012[ ) with no 
coincident WR signatures. 

- The density maps presented by [Monreal-Ibero et aT] ( [2010| 
for [S ii]) and [Monreal-Ibero et al.[ ( [2012[ for [O ii]) show that 



the gas density peak (400-500 cm ^jTsalso elongated in the 
NW-SE direction, and centred to the NW of the supernebula. 
The gas pressure implied by these densities (assuming a tem- 
perature of 10^ K) is a few xlO^ cm"^ K. From radio ob- 
servations, [Turner & Beck[ ( [2004[ ) measure densities of 3- 
4x10"^ cm""^ , meaning that the pressure deep within the < 1 pc 
size supernebula is certainly much higher. 
The Ha line near the supernebula comprises a narrow/bright 
and broad/faint component, where the broad component has 
FWHM- 100-280 kms'^ The radial velocity map (Pig.^ 
and pv diagram (Fig. [8]) show that there is a large velocity 
gradient along the same NW-SE direction traced by both the 
broad and secondary narrow Ha components (C2 and C3). 
Our new data support the conclusion of [Monreal-Ibero et al.[ 
([2010 ) that this is an ionized gas outffow, with the dynamical 
centre implying that the driving source(s) are the supernebula 
clusters (despite being so heavily embedded). 
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- |Beck et aL]p012| ) observed [S iv] 10.5 jim emission from the 
supernebula, and found a similar line profile to Ha, compris- 
ing a broad and narrow component. The broad component 
was ~20 kms"^ blueshifted and interpreted as as indicating 
a gas outflow. The peak velocity of the [S iv] agrees with the 
Brackett lines ( [Turner et al.||2003| ) and with Ha (this work; 
although the Ha broad component is broader and contributes 
less to the total line flux than [S iv]). Furthermore, spatially 
resolved observations of the radio-recombination line H53Qf 
over the same region show the same velocity gradient men- 
tioned above (Rod riguez-Rico et al.||20Q7] ). These results 
suggest that Ha, Bra, H53a and [S iv] are all tracing ion- 
ized material with approximately the same dynamics, despite 
originating at difl'erent depths within the nebula. 

Since star clusters have been forming throughout the central 
regions of NGC 5253 for last -10 Myr ( [Harris et al. '^SQCM), it 
is possible that previous generations have been responsible for 
the build-up in chemical abundances. However, the fact that the 
supernebula is located approximately in the centre of the chem- 
ically enriched region, and there is no evidence for older, suffi- 
ciently massive clusters in the vicinity strongly implies that the 
supernebula cluster(s) must be the source. 

In line with above arguments regarding chemical fluctuations 
in dwarf galaxies, we propose that in NGC 5253, the chemi- 
cally processed material expelled from star clusters in their early 
evolution remains in the hot wind phase and does not mix with 
cooler gas. However, the unusually high densities and pressures 
in the supernebula and surrounding medium have acted to stall or 
impede the winds from these embedded clusters, allowing time 
for the chemically enriched hot phase to mix with the cool/warm 
phases. Cluster winds stalled by high ambient pressures have 
been previously found in M82 ( [Silich et aLl|2007t [Smith et alT 
[2006) ). 

Are the supernebula clusters old enough to be the enrich- 
ment source? The equival ent width of Paof constrains the age 
of the clusters to < 6 Myr ( Alonso-Herrero et aL][2004[ [Martin- 



Table 5. 0/H, N/H and N/0 abundance ratios from the literature and 
this work for various regions across the starburst. The nomenclature is 
from the respective papers. These values were used to create the plots 
inFig.[T0| 



[Hernandez et al.[2005| ). That the radio and millimetre continuum 
emission are almost entirely due to thermal emission ( Beck et aL| 
[1996 ; Tu rner et al.[[2000| ) suggests that no supemovae have ex- 
ploded in the clusters, implying an upper age limit of ~4 Myr. 
Martm-Hernandez et al.[ ( [2005| ) modelled the mid-IR spectra of 
the embedded SSCs and found two solutions consistent with 
the observations: 1) a young (<4 Myr) cluster with a "non- 
standard", low, upper mass cutoff of the IMF; or 2) a cluster of 
-5-6 Myr with a "standard" high upper mass cutoff. The pres- 
ence of N-enrichment strongly argues for the influence of WR 
stars, whose appearance in a cluster peaks between ~3-5 Myr 
([Crowther 2007 ). As Martm-Hernandez et al. ( 2005 ) comments, 
one expects a supernova rate of at least 100 SN/Myr after ap- 
proximately 3 Myr, but no effects of supernovae are seen in the 
radio. We therefore suggest that the clusters must be ~4 Myr old 
- old enough to have (had) a significant WR star population, but 
young enough to have had none or very few SNe. 

That the thermal radio emission from the supernebula is op- 
tically thick classifies it as an "ultra-dense" H ii region with pres- 
sures significantly higher than typical interstellar values ( [Turner 
[et al.[[T998] ). Under normal circumstances this would imply 
that the clusters are extremely young (<1 Myr) from pres- 



sure/lifetime arguments (cf . [Kobulnicky & Johnson|1999[[Vacca 



[et al.[|200 2). However, we now know that the clusters exist in 
a much larger high density/pressure environment, helping to ex- 
plain why they remain so embedded even at ages of a few Myr 
(cf. [Smithetal.|2006D . 





arcsec(pc) 


12+log(0/H) 


12+log(N/H) 


log(N/0) 


Walsh & Roy (1989) 








Region 2 


6.68(123) 


8.22 ± 0.09 


6.87 ± 0.05 


-1.35 ± 0.05 


Region 4 


9.09(167) 


8.22 ± 0.08 


6.78 ± 0.05 


-1.43 ±0.05 


Region 5 


6.90(127) 


8.12 ±0.06 


6.80 ± 0.04 


-1.32 ±0.04 


Region 6 


7.83(144) 


8.05 ± 0.07 


6.82 ± 0.04 


-1.23 ±0.04 


Region 7 


12.70(234 ) 


7.96 ± 0.07 


6.65 ± 0.04 


-1.31 ±0.04 


Kobulnicky et al. (1997) 








HII-1 


1.50(28) 


8.12 ±0.05 


b 


-0.83 ± 0.09 


HII-2 


1.75(32) 


8.19 ±0.07 


b 


-0.86 ±0.10 


UVl 


2.70(50) 


8.16 ±0.12 


b 


-1.29 ±0.13 


Lopez-Sanchez et al. (2007) 








Knot A 


1.57(29) 


8.45 ± 0.12 


7.54 ± 0.14 


-0.91 ± 0.18 


KnotB 


0.76(14) 


8.37 ±0.13 


7.40 ±0.15 


-0.97 ±0.19 


KnotC 


2.80(52) 


8.54 ±0.11 


7.08 ±0.15 


-1.50 ±0.08 


KnotD 


7.20(132) 


8.31 ±0.07 


6.82 ± 0.07 


-1.49 ±0.10 


Monreal-Ibero et al. (2012) 








Knot 1 


0.00(0) 


8.24 ± 0.09 


7.43 ± 0.07 


-0.81 ± 0.25 


Knot 2 


2.80(51) 


8.32 ± 0.04 


6.94 ± 0.03 


-1.37 ±0.11 


Knots 


7.21(133) 


8.26 ± 0.08 


7.00 ± 0.06 


-1.26 ±0.22 


This work - 


Position 3 








Region 1 


7.4(136) 


8.40 ±0.10 


6.70 ± 0.06 


-1.70 ±0.12 


Region 2 


8.3(153) 


8.16 ±0.18 


6.61 ±0.12 


-1.55 ±0.21 


Region 3 


9.5(175) 


8.11 ±0.21 


6.68 ±0.15 


-1.43 ±0.26 


Region d_3 


9.6(177) 


8.27 ± 0.04 


6.79 ± 0.04 


-1.48 ±0.05 


This work - 


Position 4 








Region 4 


10.9(200) 


8.40 ± 0.02 


6.86 ± 0.02 


-1.53 ±0.03 


Region 5 


10.3(189) 


8.65 ± 0.34 


7.04 ± 0.20 


-1.61±0.39 


Region 6 


13.3(245) 


8.20 ±0.13 


6.58 ± 0.07 


-1.62 ±0.15 


Region d_4 


11.7(215) 


8.42 ±0.14 


6.91 ±0.11 


-1.52 ±0.18 



' Projected distance from the supernebula. 

' These authors measured the N/0 ratio directly from the intensity of [Nii]/16584 
and [O n]/l3727 since they did not detect sufficient lines to measure the N abun- 
dance itself. 



The elongated morphology of the enriched region is spatially 
coincident with the extent and orientation of the outflow, imply- 
ing that this has acted to increase the efficiency of either the mix- 
ing (to bring the new chemicals into the observable warm phase), 
and/or the subsequent diff'usion. The outflowing N-enriched ma- 
terial has an extent of ~60 pc, which could have propag ated there 
over a few Myr with difl'usion speeds of 20-60 km s"^ (Monreal- 



|Ibero et alT 2010 ). This is in good agreement with the the outflow 
speeds measured above. 

Further observations are required to explain why the densi- 
ties/pressures in this supernebula region are so high. One possi- 
bility is that the difl'use H i cloud infalling along the minor axis 
of the galaxy found by Lop ez-Sanchez et aL] ( [2012| ) may be re- 
sponsible for compressing the gas at this point. 

It is interesting that the situation in the NGC 5253 superneb- 
ula is remarkably reminiscent of that in the only other system 
with a confirmed N/H enhancement: Mrk996 (D^22 Mpc). Here 
the N/H enhancement is seen in the nuclear regi on, spatially 
coincident with a significant WN star population ( [James et"ar 
2009[). Although the elevated N/H ratios are only found in a 



broad emission-line component (FWHM>300 km s ^), the asso- 
ciated gas densities are also extremely high (10^-10^ cm"^). 



6. Conclusions and Summary 

We obtained Gemini-S/GMOS-IFU observations of four regions 
near the centre of the nearby dwarf starburst galaxy NGC 5253. 
We identify 1 1 distinct WCE-type stars through mapping out the 
strength and location of the red C iv WR bump. Some are as- 



sociated with known young clusters catalogued by [Harris et al. 
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IFU position 3: 


Region 1 


Region 2 


Re^ 


^ion 3 


Diffuse Region 




Fobs 


Fdered 


Fobs 


Fdered 


Fobs 


Fdered 


Fobs 


Fdered 


[On] ^3727+29^ 


308.30±3.40 


368.59+5.08 


308.30+3.40 


403.85+5.86 


308.30+3.40 


413.97+13.49 


308.30+3.40 


376.05+5.18 




100.00±0.89 


100.00+0.89 


100.00+0.82 


100.00+0.82 


100.00+3.10 


100.00+3.10 


100.00+0.93 


100.00+0.93 


[Oiii] /15007 


523.74±3.32 


510.84+3.29 


361.04+2.11 


347.68+2.08 


373.89+8.24 


358.81+8.06 


336.34+2.20 


327.14+2.18 


[Fern] /15270 


0.63±0.03 


0.59+0.03 


0.88+0.12 


0.79+0.10 


1.81+0.18 


1.61+0.16 


0.87+0.04 


0.81+0.04 


[Nn] ^5755 


0.30±0.03 


0.26+0.03 


0.51+0.12 


0.42+0.09 


0.77+0.11 


0.62+0.09 


0.53+0.02 


0.46+0.02 


He I ^5876 


13.82±0.13 


11.90+0.14 


14.38+0.10 


11.47+0.12 


15.82+1.15 


12.36+0.95 


13.21+0.09 


11.18+0.11 


[Oi] ^6300 


4.69±0.04 


3.88+0.05 


6.32+0.09 


4.76+0.09 


10.91+0.41 


8.00+0.39 


8.05+0.06 


6.52+0.08 


[Siii] ^6312 


2.46+0.06 


2.04+0.06 


2.22+0.09 


1.67+0.07 


2.34+1.10 


1.71+0.81 


2.11+0.04 


1.71+0.03 


Hot 


380.36±2.58 


303.62+3.77 


399.67+2.33 


284.30+3.77 


413.20+9.12 


284.88+12.69 


365.38+2.39 


284.38+3.50 


[Nil] ^6584 


19.89±0.49 


15.67+0.42 


29.88+0.29 


20.83+0.33 


33.17+1.14 


22.38+1.20 


31.11+0.22 


23.86+0.31 


He I ^6678 


4.32±0.07 


3.42+0.07 


4.17+0.15 


2.92+0.11 


5.74+0.35 


3.90+0.29 


3.90+0.06 


3.00+0.06 


[Sii] ^6716 


31.23±0.22 


24.60+0.32 


43.94+0.28 


30.64+0.43 


62.01 + 1.43 


41.83+1.96 


53.02+0.35 


40.66+0.52 


[Sii] ^6731 


23.67+0.20 


18.64+0.26 


34.43+0.23 


24.00+0.34 


45.75+1.09 


30.86+1.46 


39.19+0.26 


30.05+0.39 


dm 


0.303+0.014 


0.458+0.016 


0.501+0.052 


0.337+0.014 


erg s"^ cm"^ 


15.75+0.20 


4.67+0.010 


1.45+0.1 


12.39+0.20 


IFU position 4: 


Region 4 


Region 5 


Region 6 


Diffuse Region 




Fobs 


Fdered 


Fobs 


Fdered 


Fobs 


Fdered 


Fobs 


Fdered 


[On] A3121+29^ 


292.90+3.10 


352.87+10.08 


292.90+3.10 


335.43+5.03 


292.90+3.10 


336.22+4.77 


292.90+3.10 


357.27+5.98 


UJ3 


100.00+2.55 


100.00+2.55 


100.00+0.83 


100.00+0.83 


100.00+0.92 


100.00+0.92 


100.00+1.34 


100.00+1.34 


[Oni] /15007 


346.71+6.32 


337.81+6.29 


432.65+2.54 


424.54+2.57 


376.21+2.45 


369.03+2.45 


339.86+3.23 


330.56+3.20 


[Feni] /15270 


0.67+0.36 


0.62+0.34 


0.78+0.12 


0.74+0.11 


0.32+0.04 


0.30+0.04 


0.72+0.06 


0.67+0.06 


[Nn] /15755 


-+- 


-+- 


0.28+0.12 


0.25+0.11 


0.31+0.04 


0.28+0.04 


0.37+0.06 


0.32+0.05 


He I ^5876 


13.47+0.30 


11.53+0.36 


13.53+0.10 


12.08+0.14 


12.36+0.09 


11.01+0.12 


13.58+0.14 


11.49+0.17 


[Oi] ^6300 


9.69+0.32 


7.96+0.35 


5.68+0.09 


4.92+0.10 


2.31+0.04 


2.00+0.04 


6.82+0.08 


5.53+0.10 


[Sm] ^6312 


1.83+0.27 


1.50+0.23 


2.63+0.09 


2.28+0.08 


2.00+0.03 


1.73+0.03 


2.11+0.05 


1.71+0.05 


Ua 


362.26+6.60 


286.39+10.91 


344.89+2.03 


290.66+4.25 


337.15+2.19 


283.30+3.84 


378.06+3.59 


294.24+5.57 


[Nn] ^6584 


29.01+0.89 


22.62+1.06 


21.77+0.27 


18.16+0.34 


16.92+0.14 


14.07+0.21 


29.60+0.31 


22.70+0.46 


He I ^6678 


4.17+0.46 


3.27+0.38 


3.89+0.15 


3.25+0.13 


3.81+0.06 


3.18+0.06 


4.11+0.09 


3.17+0.09 


[Sn] ^6716 


54.72+1.05 


42.66+1.72 


31.78+0.22 


26.51+0.42 


20.76+0.14 


17.26+0.25 


50.66+0.49 


38.85+0.77 


[Sn] ^6731 


39.37+0.80 


30.70+1.25 


26.06+0.19 


21.74+0.35 


15.29+0.11 


12.72+0.18 


37.12+0.36 


28.47+0.56 




0.316+0.045 


0.230+0.018 


0.234+0.016 


0.337+0.022 


F(BJ3) X 10-15 


1.31+0.10 


4.43+0.10 


14.33+0.2 


10.75+0.5 


erg s^i cm""^ 



















^ Values assumed from i Walsh & Rovl ll989 1, their region 2. 
^ Values assumed from Walsh & Roy 1 1989 1, their region 7. 

Table 3. Integrated emission line fluxes for the six defined regions within our IFU positions 3 and 4 
H/? fluxes are given in the last rows), and de-reddenened using the c(Hy5) values shown at the bottom 
decrement, as described in the text. 



. Values are given relative to F(H/?)=100 (the 
. c(H/?) was derived from the Ha-HjS B aimer 



2004| ), and two are coincident with the location of WN-type WR 
stars detected by |Monreal-Ibero et al.| ( [2Q10} . The finding of WR 
stars spread out over 20'' (-350 pc) attests to the large area over 
which the starburst has occurred. 

We fitted multiple Gaussian profile models to the emis- 
sion lines to recover kinematic information about the ionized 
gas in the four IFU positions. Together with the results of 
Monreal-Ibero et aT] p010| ), our results paint a picture of the 
gas dynamics in NGC 5253 consisting of localised gas flows, 
as part of multiple overlapping bubbles and filaments driven 
by the star clusters throughout the starburst. We identify a 
broad (FWHM>100 kms"^) line component; near the famous 
supernebula ( Turner et al.|2003| it contributes 10-50 percent to 
the total line flux, but this falls ofl' with distance, and is not de- 
tected at all in IFU position 4 (>10'' to the south). This suggests 
that its origin is connected to the concentration of young clusters 
in the supernebula region. 

We confirm the presence of a strong velocity gradient over 
~4V5 (- 80 pc), centred on th e locati on of the supernebula, as 
seen by |Monreal-Ibero et aL] ( |2010| ), and show that both the 
broad and secondary narrow Ha components follow the same 
pattern. In line with |Monreal-Ibero et al. ( |2010 |), we interpret 
this as an accelerating ionized gas outflow from the supernebula 
clusters. 

We extracted integrated spectra and measured ionized-phase 
abundances of O, N and He in a number of regions at distances 
7V4-13V3 from the supernebula. By combining these with 0/H, 
N/H and N/0 measurements of nearer regions found in the liter- 
ature, we found that, within the central 250 pc, the 0/H and N/H 



distributions are flat and consistent with a metallicity of -0.2 Z©. 
It is only within the central 50 pc that the nitrogen abundance is 
strongly enhanced by log(N/H)~0.5 dex. 



This enhancement is what puts NGC 5253 in the rare class 
of systems with significant chemical variations. That the fa- 
mous supernebula is located approximately in the centre of the 
N-enriched region, and there is no evidence for older, sufficiently 
massive clusters in the vicinity, strongly implies that the su- 
pernebula cluster(s) must be the source of the extra N. The pres- 
ence of this enrichment strongly argues for the influence of WR 
stars and thus an age for the supernebula clusters of ~3-5 Myr; 
the lack of evidence for supemovae argues for ages <4 Myr 
dBeck et al.|1996t|Tumer et al.|2000| ). 

We propose that the unusually high densities and pressures 
in the supernebula and surrounding medium have acted to im- 
pede the winds from the supernebula clusters, allowing time 
for the chemically enriched hot wind material to mix with 
the cool/warm phases (and therefore become visible in the 
optical/near-IR). Under normal conditions, this enriched wind 
material remains in the ho t wind phase and does no t mix with 
cooler gas (as described by |Kobulnicky & Skillman|1997^ ), thus 
explaining why the rest of the galaxy (and most other dwarf 
galaxies) is chemically uniform in the warm ionized gas phase. 
We argue that the outflow has acted to spread the injected N-rich 
gas from the WR winds and helped to increase the efficiency of 
the hot-cool phase mixing to allow it to be seen in the optical. 
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Quantity 


Region 1 


Region 2 


Region 3 


Pos.3 Diffuse 


Region 4 


Region 5 


Region 6 


Pos.4 Diffuse 


T.([Nii]) 
T.([Om])^ 


96±17 
10555±500 
10194±906 


153±17 
10564±1035 
11492±1639 


58±42 
11775+1288 
11486+2102 


56+12 
11187+219 
10651+382 


35+30 
9700+0 " 
10000+0 


208+48 
8048+883 
7471+1323 


57+20 
11378+758 
10639+927 


51+20 
9681+769 
8850+815 


0+/i^+ X 10^ 
X 10^ 

O/HxlO^ 
12+log(0/H) 


1.06±0.19 
1.45±0.34 
2.51±0.56 
8.40±0.10 


0.94±0.34 
0.54±0.27 
1.46±0.60 
8.16±0.18 


0.56+0.25 
0.63+0.38 
1.30+0.64 
8.11+0.21 


0.88+0.07 
0.92+0.10 
1.86+0.17 
8.27+0.04 


1.33+0.07 
1.16+0.04 
2.50+0.10 
8.40+0.02 


2.35+1.48 
1.84+1.93 
4.46+3.51 
8.65+0.34 


0.65+0.16 
0.94+0.27 
1.58+0.47 
8.20+0.13 


1.23+0.40 
1.48+0.52 
2.65+0.89 
8.42+0.14 


X 10^ 
N/HxlO^ 
12+log(N/H) 
log(N/0) 


2.09±0.29 
5.01±0.75 
6.70±0.06 
-1.70±0.12 


2.66±0.76 
4.07±1.10 
6.61±0.12 
-1.55±0.21 


2.40+0.66 
4.81 + 1.69 
6.68+0.15 
-1.43+0.26 


3.09+0.16 
6.22+0.51 
6.79+0.04 
-1.48+0.05 


3.88+0.14 
7.32+0.33 
6.86+0.02 
-1.53+0.03 


4.65+2.16 
10.90+5.03 

7.04+0.20 
-1.61+0.39 


1.68+0.28 
3.83+0.65 
6.58+0.07 
-1.62+0.15 


3.66+0.84 
8.06+1.99 
6.91+0.11 
-1.52+0.18 


X 10^ 
X 10^ 

S/HxlO^ 

12+log(S/H) 

log(S/0) 


0.75±0.08 
3.39±1.02 
4.36±1.27 
6.64±0.13 
-1.76±0.16 


0.94±0.23 
1.56±0.98 
2.32±1.11 
6.37±0.21 
-1.80±0.27 


1.03+0.26 
0.77+0.75 
1.89+1.07 
6.28+0.24 
-1.84+0.32 


1.21+0.06 
2.94+0.46 
4.34+0.55 
6.64+0.05 
-1.63+0.07 


1.62+0.04 
3.39+0.51 
5.23+0.54 
6.72+0.04 
-1.68+0.05 


1.75+0.61 
5.89+9.62 
8.48+11.30 
6.93+0.57 
-1.72+0.67 


0.47+0.07 
2.39+0.86 
3.03+0.99 
6.48+0.14 
-1.72+0.19 


1.37+0.25 
5.88+2.63 
6.95+2.86 
6.84+0.18 
-1.58+0.23 


He/HxlO^ 


8.13±0.09 


7.99+0.05 


8.99+0.55 


7.91+0.07 


8.19+0.29 


8.19+0.29 


7.92+0.07 


7.92+0.11 


Weq(H^) (A) 

L(WR bump) 
(xlO^^ ergs-i) 

Equivalent no. of 
WCE'^ stars 


-96.55±5 
10.65±1.17 

~1 


-60.12+4 


-14.06+2 
9.05+1.17 

~1 


-55.43+3 


-33.52+5 
4.91 + 1.17 

~1 


-105.35+10 
8.49+1.17 

~1 


-99.88+10 


-42.41+5 



^ [Nii]/15755 was undetected in this region, hence we adopted the average Tg from the surrounding regions of ~9700 K. 
^ Derived from Perez-Montero relationship. See text for details. 

^ Early-type WC stars (WCE) are expected to dominate at the metallicity of NGC 5253. To estimate the number of WC stars in each region, we used the relation 
given by Lopez-Sanchez & Esteban i 2010 i for the Civ/15808 luminosity of a WCE star [i.e. L(WR bump)] at the metallicity of NGC 5253. The luminosities we 
measure are consistent with the presence of 1 WCE star in each source. 

Table 4. Abundance measurements for the six defined regions within IFU positions 3 and 4. Also given are the equivalent widths of HjS (Weq(HyS), 
and the integrated red WR bump luminosities that provide an estimate of the equivalent number of WCE stars. 
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